Condensation of helium in nanoscopic alkali wedges at zero temperature 
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We present a complete calculation of the structure of liquid ''He confined to a concave nanoscopic 
wedge, as a function of the opening angle of the walls. This is achieved within a finite-range density 
functional formalism. The results here presented, restricted to alkali metal substrates, illustrate the 
change in meniscus shape from rather broad to narrow wedges on weak and strong alkali adsorbers, 
and relate this change to the wetting behavior of helium on the corresponding planar substrate. As 
the wedge angle is varied, we find a sequence of stable states that, in the case of cesium, undergo 
one filling and one emptying transition at large and small openings, respectively. A computationally 
unambiguous criterion to determine the contact angle of ''He on cesium is also proposed. 

PACS numbers: 68.08.Bc, 67.70.+n, 68.03.-g, 67.90. +z 
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I. INTRODUCTION 

The current possibility of tailoring the production of 
adsorbing surfaces in the nanoscopic scale, raises a series 
of questions on the role of inhomogeneities in the growth 
of liquid films on nanopatterned planar substrates. In 
particular, the wedge geometry is a simple one that per- 
mitted to carry on many theoretical anticipations of the 
wetting behavior of fluids in such confinement (see, e.g., 
Refs. and cited therein), and it has been recently 

argued^, that models of fluids adsorbed in planar wedges 
can be regarded as a generic class of structured inhomo- 
geneous fluids. 

Since the early calculation of physisorption of noble 
gases in an oblique corner f4 most predictions of the shape 
changes undergone by the liquid meniscus are based on 
the examination of the equilibrium liquid-vapor inter- 
face of a macroscopic sample at constant bulk density 
POi whose sharp free surface is described by a function 
determined by minimization of the surface free energy. 
The latter contains the energy cost for building a curved 
liquid- vapor interface with surface tension aiy, in ad- 
dition to the interfacial contributions governed by the 
surface energies a si and dav, with l,s,v denoting liquid, 
solid and vapor/ vacuum, respectively. These three sur- 
face tensions determine the contact angle 9 as given by 
the classical Young's relation These calculations con- 
firm the thermodynamic theory of wetting in a wedged 
that shows, on the ground of purely macroscopic argu- 
ments, that a filling transition characterized by a per- 
fectly planar liquid- vapor interface takes place, at liquid- 
vapor coexistence, when the opening angle 2lp of the walls 
is related to the contact angle for a planar surface as 
tf — 71/2 — 6. Moreover, the above theories classify the 



meniscus shapes for nonwetting liquids according to the 
following sequence 
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Lp, spherical bridge; 
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(f, spherical convex meniscus; 



• 6 — ^ — Lp, planar meniscus; 

• 6 < ^ — ip, spherical concave meniscus. 

While classical, macroscopic approaches provide a solid 
frame for the expected behavior of fluids in wedges, one 
may wonder if these predicted phenomena are robust 
at the nanoscale. A study of interfacial phenomena in 
power-law wedges permits to trace the route from wet- 
ting to capillary condensation^ and suggests the possibil- 
ity of adapting the adsorptive abilities of solid substrates 
by shaping its surface inhomogeneities at the nanoscopic 
level. Recent measurements of the growth of Ar films 
on an array of microscopic linear wedges^ confirm this 
conjecture and demonstrate a clear crossover between a 
planarlike and a geometry-dependent behavior. On the 
theoretical side, it has been remarked^ that the full one- 
body density p(r) for the confined particles should be 
obtained either from simulations or from density func- 
tional (DF) theory. With the latter instrument, we have 
recently investigated the growth of nanoscopic ''He clus- 
ters on planar alkali surfaces at zero temperature (T), 
employing a finite range DF.l'^ The analysis of the ener- 
getics and shape systematics with increasing number of 
atoms N shows that helium drops on wettable alkali sur- 
faces reach a maximum vertical height equal to the thick- 
ness of the prewetting film, with the chemical potential at 
the prewetting jump in the zero temperature adsorption 
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isotherm. These results iUustrate that DF-based calcu- 
lations can account for the interplay between geometry 
and thermodynamics in nanoscopic systems, sufficiently 
large to remain beyond the reach of present microscopic 
many-body theories. On these theoretical grounds, in 
this work we carry out a systematic investigation of the 
spatial structure and shape evolution of nanoscopic he- 
lium samples in linear alkali wedges, as described in the 
next Sections. In particular, in Sec. ^we shortly review 
the DF formalism employed, and in Sec. 11111 we present 
our results for Cs and Na wedges, which are summarized 
in Sec. HVl 



II. METHOD 

Our geometry is translationally invariant along one 
spatial direction (y-axis). The equilibrium "^He density 
profile p{x, z) will thus depend on the {x, z) coordinates 
only. We compute the density profile that minimizes the 
zero temperature grandpotential ~ E—fi N, with /i the 
chemical potential of the helium atoms and with TV the 
grandcanonical ensemble average of the particle number 
operator, that in a wedge extending a macroscopic length 
L along the y-direction, takes the form 



dx dz p{x, z) 



(1) 



The variation 6fl/Sp gives rise to the integrodifferential 
Euler-Lagrange equation 
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FIG. 1: (Color online) Chemical potential of a *He pancake 
on planar Cs as a function of linear density. The chemical 
potential of the n — 280 pancake is indicated by an 

horizontal dashed line. The inset shows, for the n = 140 
density profile, the equidensity lines from p = 0.9po to O.lpo 
in O.lpo steps. 



Equation ^ is discretized with 13-point formulas 
for the partial derivatives and solved by an imagi- 
nary time method^ employing a Fast Fourier Trans- 
form algorithmic to obtain V{p) from the atom density 
p{x, z)pi& We have verified that the numerical outcome is 
stable against the increase of the number of mesh points 
and the order of the discretization formulas. 
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+ V{p) + Vs{x,z) 



\J p{x,z) = p.^Jp{x,z) 
(2) 

Here V{p) is the effective potential arising from func- 
tional differentiation of the potential energy density per 
unit length, and Vs{x,z) is the wedge potential. In this 
work, we have selected the full Orsay-Trento (OT) DF,^^ 
which is known to provide a quite accurate description 
of inhomogeneous structures of "'He at T = 0. We can 
safely neglect gravitational effects in our calculations, 
since the capillary length a — [2aiy / (gm^pa)]^/'^ , with 
po = 0.0218 A"'^ the saturation density of ^He and g the 
acceleration of gravity, is much larger than any linear 
dimension of the "^He samples considered in the present 
work. 

To our knowledge, the best available adsorbing po- 
tential for helium on alkali planar surfaces is the 
Chizmeshya-Cole-Zaremba (CCZ) potential.^ As a com- 
promise between rigor and physical insig ht^ in the 
present work we have chosen to approximate Vs{x, z) by 
the summation of two planar CCZ potentials at an angle 
2Lp. Double-counting effects, even when sizeable, are only 
important for fairly large values of 2(p and in the vicinity 
of the wedge apex.i Thus, one should not expect qualita- 
tive differences in the density patterns for the fairly large 
linear densities considered here. 



III. RESULTS 



A. Cesium wedges 

We address first the most interesting case of cesium, 
which is not wetted by *He at T = 0. For reference, we 
have computed a sequence of ^He pancake-shaped sys- 
tems on a ffat Cs surface, namely, translationally invari- 
ant systems in the ^/-direction and characterized by the 
linear density n = N/L. Fig. ^ shows the chemical po- 
tential of the ''He atoms as a function of n, that displays 
a neat tendency towards saturation at the bulk figure 
/ig — —7.15 K. The equidensity lines for the density pro- 
file corresponding to a linear density n = 140 A"'^ are 
shown in the inset of this figure. The profi les are quali- 
tatively similar to those shown in Ref. [13 in the case of 
the deposited drops, and confirm the nanoscopic spatial 
scale of the systems. 

It should be kept in mind that in a typical adsorption 
experiment, one changes the chemical potential by vary- 
ing the pressure in the chamber towards saturated vapor 
pressure (SVP) conditions. Very close to SVP (but below 
it) the thermodynamic equilibrium state of liquid ^He on 
Cs is a thin, microscopic planar film partially wetting 
the surface. Nonetheless macroscopic, metastable ''He 
droplets can be realized if a finite amount of helium is 
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FIG. 2: (Color online) Density profiles of helium at n = 280 
in a Cs corner represented by thin straight lines, for 
angular opening 2ip = 15, 45, 52, 75, 110, and 120°. The 
distance between major tick marks on the z axis is 50 A. The 
center of the 15° configuration is about 400 A above the apex 
of the wedge. Equidensity lines are drawn as in Fig. Q 



deposited on the Cs surface, as shown for instance in 
the experimental results of Ref. The contact angle of 
such droplets with the surface plane is a meaningful quan- 
tity to be measured/calculated and which can be related 
to thermodynamic -equilibrium- properties. We shall see 
in the following that similar non-equilibrium configura- 
tions are found in the wedge geometry, which might be 
experimentally realized under particular conditions. 

When helium at given n is placed in a Cs wedge and 
the opening angle 2ip is varied, we generate the shape 
sequence illustrated in Fig. El at a linear density of 280 
A"-'^. This evolution illustrates the modifications of the 
wetting behavior of helium on Cs due to the confine- 
ment exerted by the corner. A convex helium pancake 
like the one shown in the inset of Fig. ^ initially ad- 
sorbed on a flat Cs surface, when subjected to the for- 
mation of a wedge and decrease of the opening angle 
experiences a sequence of changes whose highlights are: 
a) a filling transition where the curvature vanishes, for 
an opening near 110°, b) formation of a concave menis- 
cus that persists within a sizable range of angles; c) an 



-7.05 



- \ ^ 

\ ^ ^ 

- ^ \ 


n= 70 a' 

n= 140 a' 

n = 280 A"' 

Cs 




- y^ 

' ^^^ — " / 

" 1 , , 1 , , 1 , , 1 , , 1 , , 1 , , 1 , , 1 



J I 1 I I I I I I I I I I I I I I I I I I I 

15 30 45 60 75 90 105 120 
2(p (deg) 



FIG. 3: (Color online) Chemical potential of helium atoms in 
a Cs wedge as a function of angular opening 2<^ (degrees) for 
linear densities 70, 140 and 280 A~^. Wedge configurations 
to the left of the thin vertical segment at about 43° have 
higher energies than a pancake, with the same linear density, 
adsorbed on a planar surface. 



'emptying' transition where a flat meniscus appears near 
52°, with restoration of the convex shape for smaller an- 
gles; d) formation of a bridge that will be expelled from 
the wedge for a sufflciently small angle. We have verified 
that the bridge profiles are perfectly circular, as predicted 
by the classical theories A^iSii Note, however, that here 
the term 'transition' does not address a thermodynamic 
phase transition, but a crossover between two different, 
well defined regimes characterized by geometry and en- 
ergetics as analyzed in what follows. 

The above sequence can be understood by examining 
the energetics of helium atoms at several linear densities, 
as functions of the wedge opening. For this sake, in Fig. 
Olwe plot the chemical potential for n = 70, 140 and 280 
A"-'^. Several interesting features become evident. First, 
the three curves cross at the bulk value /iq for two open- 
ing angles, namely the filling one 2(pp = 110° and the 
emptying one 2(pE = 52°. Within this interval, the he- 
lium chemical potential lies below the thermodynamic 
limit fiQ, and the systems are thermodynamically stable, 
since dfi/dn is positive. By contrast, above the filling 
angle Lpp and below the emptying angle ipE, fJ- lies above 
the bulk value with negative dfi/dn; the corresponding 
helium samples are then metastable. These characteris- 
tics can be further analyzed looking at the grandpotential 
per unit length fl/ L = E/L — iin, where E/L is the en- 
ergy per unit length, displayed in Fig. ^ for the same 
linear densities. Note that in the present context, il/L 
represents the surface-plus-line energy difference per unit 
length between a wedge filled with liquid and an empty 
wedgeii At (pp, ^/L is independent of n and vanishes 
for these large, however nanoscopic, amounts of fluid. 
Thus, the conditions for a filling transition, previously 
derived for macroscopic systems,- are fulfilled for these 
helium samples. This suggests an effective way of deter- 
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FIG. 4: (Color online) Same as Fig. |21for the grandpotential 
per unit length. 

mining the contact angle B of ^He on Cs by means of 
DF calculations; we recall that previous DF-based cal- 
culations of this angle give values between 31° (Ref. [l7|) 
and 36°fi& in good agreement with experiments in Refs. 
ma Other DF estimates provide somewhat smaller 
values f-- whereas a different experimental determination 
of the contact angle^^ yields Q — 48° ± 1°. Using the 
value of i-pp ^ 110° where the curves cross in Fig. 0]we 
obtain here a contact angle 9 = n /2 — tfp = 35°. 

Figure 21 also shows, in correspondence with Fig. 
13 that a second curve crossing, where 17 /L is n- 
indcpendent, occurs at the aperture 2ipE ~ 53° i.e., at 
the threshold for the formation of the bridge, which we 
have found to be clearly visible only below ~ 30°. Except 
in a very narrow angular region near (f e , whose existence 
is likely attributable to finite size effects^^ Vt/L is neg- 
ative in the interval [(PEt^Pf], indicating that condensa- 
tion might occur either in a continuous way or through 
discontinuous jumps in the linear coverage. In fact, the 
latter is the case here, as shown below in relation to Fig. 

n 

Within the interval [lpet^f], the gain in energy origi- 
nates in the disappearance of the solid-vacuum interface 
at the meniscus contact line, that leaves only a liquid- 
vacuum interface at the free surface, as seen in Fig. |21 
This is reversed for p < LpE since the system shifts to a 
nonwetting regime. In these wedges, the energy balance 
favors the presence of vacuum, rather than liquid, at the 
corner, making room to the formation of the inner convex 
meniscus that 'dries' the corner and forms the bridge, at 
the price of increasing the total energy. 

The existence of two stable solutions in the [pet'^f] 
interval for a given coverage as seen in Fig. and in 
flat meniscus configurations at = /io, is a genuine 
nanoscopic effect, as it docs not appear for macroscopic 
samples, where only a filling transition is predicted for 
liquid in a wedge made of partially wettable surfaces.^ 
As a matter of fact, all equilibrium macroscopic config- 
urations, irrespective of the opening, might be thought 
of as corresponding to the same chemical potential, that 
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FIG. 5: (Color online) Chemical potential of helium atoms 
in a Cs wedge as a function of linear densities. Solid line: 
2ip = 49°; dashed line: 2ip = 53°. 



of bulk liquid, with the fluid density in the wedge tak- 
ing its value at bulk liquid-vacuum coexistence at zero 
temperature. 

To illustrate the mechanism of condensation, in Fig. [51 
we plot two selected adsorption isotherms ^{n), for an- 
gles 2ip = 49 and 53°. The respective n = ordinates 
-not shown in the current scale- correspond to the bind- 
ing energy of a single helium atom to the wedge at each 
aperture, computed separately. The solid line {2p — 49°) 
represents ^{n) for an opening slightly below LpE'- the 
^He configurations in such wedge are always metastable 
with respect to the bulk, since /i(n) approaches /iq from 
above. Thus, no condensation can occur in the wedge at 
such an aperture. By contrast, when the angle is slightly 
larger that tpE (dashed line, 2tp = 53°), condensation 
takes place with a jump in the linear density from n — Q 
to a finite value. The associated 'pre-emptying' transi- 
tions occur at a value of the chemical potential which can 
be found by using the usual equal-area Maxwell construc- 
tion (or, equivalently, by identifying the crossing between 
the ii{n) and E/N curvesiiS) From this construction we 
find that at 2p — 53° a jump from n = to n ~ 270 
occurs ai fi — — 7.155 K, i.e. just below SVP. This in- 
dicates that a crossover between a filled and an empty 
wedge is only meaningful for linear densities near to and 
above - 270 A-\ i.e., more dilute helium vapors can- 
not condense in the wedge at angles equal to and slightly 
above 2pE- A similar behavior (i.e. finite jumps in cover- 
age below SVP) is encountered for all angles in the range 
^^PE^2^pF\', however, for the largest openings in this in- 
terval, condensation occurs at rather low coverages, say 
n « 5-10 A-i. 

The occurrence of an 'emptying' regime at lower open- 
ing angle, as suggested by our results, is surprising at 
first sight. No 'emptying' transition is predicted from 
macroscopic thermodynamic arguments, and this behav- 
ior seems also to contradict the expectation that a narrow 
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wedge should be filled by capillary condensation (CC), 
below SVP, in the same way as small pores/slits are 
known, from theoretical calculations, to undergo a CC 
transition even for weakly adsorbing surfacesi2i24 We 
believe that the wedge geometry makes here an impor- 
tant difference with respect to pore/slit case -even in the 
narrowest wedge, the distance between the walls can be 
arbitrarily large sufficiently far from the apex. 

The analysis of the energy per particle for fixed n, as 
a function of wedge opening, indicates that for the linear 
densities under consideration, below 2(/3 « 43° helium 
bridges are metastable (i.e. have higher energies) with 
respect to helium pancakes deposited on planar Cs sur- 
faces. This limit is shown in Fig. Elby a vertical segment. 
Similarly, wedges corresponding to openings 2(p > 120° 
are also metastable with respect to pancakes with the 
same number of atoms. It should be remarked that these 
metastable configurations appear for finite average num- 
bers of "^He atoms (per unit length) within a grandcanon- 
ical description; as a reference, we recall that a droplet 
with finite atom number N on Cs is metastable, not only 
with respect to bulk, but also with respect to a thin 
film at finite coverage N/A. In view of these results, we 
are confident that in an hypothetical experimental setup 
where finite, nanoscopic amount of He are deposited in 
Cs alkali wedges, our prediction that narrow wegdes are 
either empty or hosting 'bridge' configurations could in 
principle be checked. 



B. Sodium wedges 

We have carried out similar calculations for helium 
on sodium, analyzing the evolution of helium initially 
placed in the unstable regime at low coverages -i.e., be- 
low the prewetting transition-. Fig. [^displays the chem- 
ical potential of ^He pancakes on a flat Na surface. For 
this strong adsorber, the corresponding asymptotic limit 
is —8.75 K, the value of the chemical potential at the 
prewetting jumpiiSi 
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FIG. 7: (Color online) Chemical potential of helium atoms in 
a Na wedge as a function of angular opening 2ip (degrees) for 
a linear densities n = 10 and 15 A~^. For n = 15 A~^, wedge 
configurations within the thin vertical segments are stable 
with respect to the infinite pancake configuration. For n = 10 
A~^, wedge configurations to the right of the corresponding 
thin vertical segment are stable with respect to the infinite 
pancake configuration. 
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FIG. 8: (Color online) Same as Fig. Qfor the grandpotential 
per unit length. 
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FIG. 6: (Color online) Same as Fig. 0for *He on planar Na. 
In this case, the inset corresponds to n = 10 A~^. 



Figs. [7| and |S1 display, respectively, the chemical poten- 
tial and grandpotential per unit length of a ''He sample 
with n = 15 A~^^ and shows that the filling transition is 
obviously absent due to the fact that sodium is wetted 
by ''He at T = 0. In these fi gures, the thresholds for the 
stability of pancake structures are shown by the short 
vertical segments. 

Finally, typical density patterns are displayed in Fig. 
Ofor the linear density n = 15 and several wedge 
openings. The panels show the spread of the wetting liq- 
uid, that retains a rather deep concave meniscus. This 
figure shows that helium originally deposited on a flat, 
strongly wettable adsorber like Na, just follows the fold- 
ing of the surface with angle 2ip and forms a deep menis- 
cus. For sufficiently small angles the pattern evolves into 
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FIG. 9: (Color online) Density profiles of helium at n = 15 
in a Na corner represented by thin straight lines, for 
angular opening 2ip — 60, 90, 105, and 120°. The distance 
between major tick marks on the z axis is 20 A. Equidensity 
lines as in Fig. |21 

a nanoscopic bridge, similar to that shown in the first 
panel of Fig. |21 Such configurations have, however, 
higher energies than the corresponding pancake struc- 
tures with the same linear density, deposited on one of 
the walls, and thus should not occur in sodium or any 
other wettable wedge. 

IV. SUMMARY 

We have presented the first detailed calculations of the 
structure and energetics of nanoscopic samples of ^He 
in a linear wedge with two identical alkali walls. Our 
results indicate that for a Cs substrate, a sequence of 
states that include one filling transition and a threshold 
for -metastable- configurations where the wedge is void 
of ''He occurs by varying the wedge opening at zero tem- 
perature, with condensation at negative grandpotential 



(accompanied by pre-filling jumps in the linear density) 
as an intermediate regime. When the material is near 
saturation, the angle for the filling transition is stable 
against changes in the linear density and can be identi- 
fied with the contact angle on a planar substrate. It is 
important to bear in mind that the experimental realiza- 
tion of such system, which could be done in practice by 
depositing the chosen alkali on a silicon surface previously 
patterned with microscale wedges is very challenging 
due to the reactive character of alkalis. 

The present wedge geometry can be also used as a 
first step to model rough alkali surfaces. Quite surpris- 
ingly, our calculations indicate that narrow Cs wedges (of 
nanoscopic size) are not filled by finite amounts of he- 
lium. This circumstance might be relevant in the study 
of the wetting behavior of "^He on rough Cs films, and in 
the modelling of helium nanopuddles formed on strong- 
pinning surfaces!^ 

We have also carried out similar calculations for he- 
lium on sodium which indicate that a helium pancake 
originally deposited on a flat, strongly wettable adsor- 
ber like Na, basically follows the folding of the surface 
with angle 2ip and forms a deep meniscus. Lastly, we 
have shown that a theoretical search of the vanishing of 
the grandpotential as a function of wedge opening, along 
with liquid- vacuum coexistence at fiQ and independent on 
the linear coverage n, is feasible within a DF description 
and unambiguously determines, from the theory point of 
view, the contact angle of a nanoscopic helium sample on 
a nonwettable substrate. 
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